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Using data collected with the DO detector at the Fermilab Tevatron Collider, corresponding to 
5.3 fb~^ of integrated luminosity, we search for violation of Lorentz invariance by examining the ti 
production cross section in lepton+jets final states. We quantify this violation using the standard- 
model extension framework, which predicts a dependence of the tt production cross section on 
sidereal time as the orientation of the detector changes with the rotation of the Earth. Within this 
framework, we measure components of the matrices (cQ)fi^33 and (€(7)^^33 containing coefficients 
used to parametrize violation of Lorentz invariance in the top quark sector. Within uncertainties, 
these coefficients are found to be consistent with zero. 

PACS numbers: ll.30.Cp, 13.85.Qks,. 14.65. Ha 



We investigate the possibility of Lorentz invariance vi- 
olation (LIV) in the top quark (t) sector, using data col- 
lected with the DO detector at the Fermilab Tevatron pp 
Collider corresponding to 5.3 fb~^ of integrated luminos- 
ity collected between August 2002 and June 2009. We 
examine events in which a ti pair is produced and decays 
into a final state including two light quarks {q,q'), two 
b quarks {b,b), and a lepton- neutrino pair {1, 1^1) via the 
mode ti -^ W^bW'b — > tvibqq'b, where £ = e,^i. The 
standard-model extension (SME) framework [1] provides 
an effective field theoretical treatment for violation of 
Lorentz and CPT symmetry in particle interactions by 
introducing Lorentz-violating terms to the Lagrangian 
density of the standard model (SM). As yet, there are 
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no quantitative limits on violations of CPT or Lorentz 
invariance in the top quark sector [2|. This parameter 
space is accessible only at high-energy particle colliders. 
Because top quarks decay before hadronizing, this study 
also offers the possibility of extending such investigation 
to what are essentially free quarks. 

Strong limits have been set on the magnitude of LIV 
in gravitational and electromagnetic interactions, as well 
as in many particle sectors. Most constraints in the par- 
ticle sectors are for matter involving quarks of the first 
generation. There are also sensitive limits on SME co- 
efficients for the second generation, but only a few for 
the third generation y. The latter include limits for 
the b quark through B-meson oscillations 3], and for r 
neutrinos from neutrino oscillations [J|. There is also a 
constraint for r leptons deduced from theoretical grounds 
using astrophysical observations |5[ . Constraints on LIV 
have also been predicted for the SM Higgs sector, as de- 
rived from radiative corrections [Q]. Many of the limits 
on LIV coefficients in the quark, lepton, and gauge bo- 
son sectors are ^10~^. However, no constraints have yet 
been placed on LIV in the top quark sector. Because the 



SME represents a general phenomenological formalism, 
the LIV terms of the SME are not constrained to cou- 
ple with the same strength to all particle species. We 
therefore consider separately only those SME terms that 
affect the top quark fields in ti events. 

While it has been shown that CPT violation implies vi- 
olation of Lorentz invariance [7| , the contributions from 
CPT violating terms in the SME to the matrix element 
for ti production and decay are suppressed. However, 
contributions from other Lorentz-violating terms can be 
significant [S]- At leading order in LIV coefficients, the 
matrix element describing the production and decay of 
a tt pair involves coefficients of the form c^^, where fi 
and V refer to space-time indices. Although at lead- 
ing order CPT-odd SME terms describing LIV in the 
top quark sector are not observable in ti production or 
decay, this analysis is sensitive to several components 
of the CPT-even {cq)^i^ab and {cu)^^ab terms, where 
A,B = 3, 3 refer to the third quark generation. The 
(cQ )/ii/33 are the SME coefficients coupling to the left- 
handed components of the third generation quark fields, 
and {cij)^u^z are the SME coefficients coupling to the 
right-handed singlet top quark field. For brevity, we 
drop the generation subscripts since we are restricting 
the analysis to the terms that couple to the top quark 
fields. To compare our results with SME studies in other 
particle sectors [2|, we also examine the linear combina- 
tions 



dfj.u — {cQ)fj.u — {cu)tJ.u- 



(1) 



The matrix element for leading-order ti production and 
decay, including leading-order contributions from SME 
terms, can be written as [Sl 



\M 



SME 



= PFF+{SP)FF + P{SF)F + PF{SF). (2) 



The P terms are functions of the parton momenta at the 
ti production vertex, while the F terms involve parton 
momenta at the decay vertices. The PFF term corre- 
sponds to the usual SM component, while the (5-terms 
reflect the dependence on SME coefficients. This expres- 
sion summarizes how the SME modifies the matrix ele- 
ment for ti production and decay at leading order. 

The i5-terms contain contractions of c^j^ coefficients 
with tensors that are functions of the four-momenta of 
the particles in ti production and decay. Due to the 
V — A structure of the weak current, the right-handed 
coefficients, (cc/)^,y, couple only to the production {6P) 
terms, while the left-handed coefficients, {cq)^v^ couple 
to both production and decay {5F) terms. The matrices 
of c^u coefficients are symmetric and traceless. Within 
the SME, these coefficients are defined by convention in 
the canonical Sun-centered reference frame [2,] . 



The kinematic component of the (5-terms of Eq. ([2]) can 
be evaluated in any coordinate system. A convenient ref- 
erence frame is that of a coordinate system fixed to the 
measuring apparatus, and we therefore choose to evalu- 
ate such contractions in the DO coordinate system. In 
this system, the momenta entering the calculation of Eq. 
([2]) are just the momenta of the particles measured in the 
detector, and, to calculate the matrix element, the coeffi- 
cients {cu)fj,u and {cq)^i, must therefore be transformed 
from the Sun's reference system to the DO coordinate 
system. 

Since the Earth is rotating about its axis, the trans- 
formation of the coefficients {cu)ij,i, and {cq)^i, from the 
Sun-centered frame to the laboratory frame introduces a 
time dependence. The relevant time scale is the sidereal 
day, which has a period of 23 hr 56 min 4.1 s (86,164.1 s). 
If any of the coefficients {cu)^i, or {cq)^i, are non-zero in 
the Sun-centered frame, they can be detected through a 
periodic oscillation in the number of ti events observed in 
the Earth-based detector as a function of sidereal time. 

The data used for this analysis correspond to 5.3 fb~^ 
of integrated luminosity collected with the DO detector. 
The DO detector [9| consists of several subdetectors de- 
signed for identification and reconstruction of the prod- 
ucts oipp collisions. A silicon microstrip tracker and cen- 
tral fiber tracker surround the interaction region for pseu- 
dorapidities jryl < 3 and \'q\ < 2.5, respectively (where 
1] = — ln[tan(0/2)] is measured relative to the center of 
the detector, and 9 is the polar angle with respect to 
the proton beam direction) . These elements of the cen- 
tral tracking system are located within a 2 T supercon- 
ducting solenoidal magnet, providing measurements for 
reconstructing event vertices and paths of charged parti- 
cles. Particle energies are measured using a liquid argon 
and uranium calorimeter. Outside of the calorimetry, 
trajectories of muons are measured using three layers of 
tracking detectors and scintillation trigger counters, with 
1.8 T iron toroidal magnets between the first two layers. 
Plastic scintillator arrays in front of the end-calorimeter 
cryostats provide measurements of luminosity. 

We employ the same event selection as described in 
greater detail in Ref. [10| . Briefly, events are collected us- 
ing a suite of triggers selecting events with a single lepton 
(e or /i) or a single lepton plus a jet. Candidate ti events 
in the lepton-|-jets channels are then selected by requiring 
the presence of one isolated electron (or muon) candidate 
with transverse momentum px > 20 GeV and pseudo- 
rapidity |?7| < 1.1 (2), and an imbalance in transverse 
energy of J^j, > 20 GeV (25 GeV). Events are divided 
into bins of jet multiplicity, and all jets are required to 
be reconstructed with pT > 20 GeV and |?7| < 2.5, with a 
leading jet of pr > 40 GeV. One of the jets is required to 
be tagged as a &-jet candidate through a neural-network- 
based (NN) algorithm llj. The time of production of 
each ti event is recorded with the event data, with an 
average accuracy of approximately ± 30 s. To follow the 



conventions utilized in other SME studies [2| , we shift the 
origin of the time coordinate to correspond to the vernal 
equinox of the year 2000. 

The SME predicts time dependent effects on the tt 
cross section of the form 



a{t) 



[1 + /sME(i)] , 



(3) 



where tTavo is the observed (time averaged) cross section 
for ti -J- W+bW-b -j- Ivfbqq'h, in £+jets final states. To 
arrive at Eq. ([3|), we compare the contribution from the 
SME terms in Eq. ^ to the SM expectation by consid- 
ering the ratio of |A^|gME ^^ the SM component PFF. 
The SME contributions in this ratio are collected into the 
function 



/sME(i) = [{cQ)t,u + {cu)^u] Kit)R^t)A 



aP 



aP 



icQ),,K{t)R-'pit)A''J 



(4) 



Eq. ((4]) is a product of the matrices of time-independent 
coefficients (cq)^i, and (cu)^iv, four-by-four matrices of 
terms that depend on the event production {Ap) and 
decay {Ap) kinematics in the DO frame, and a rotation 
matrix i?^ (i) that transforms Ap and Ap from the DO 
frame to the Sun-centered frame. 

The Ap and Ap matrices are evaluated using tt 
Monte Carlo events generated with pythia [l^j. Events 
that pass detector acceptance, trigger, event reconstruc- 
tion, and analysis selections (modeled by a full simulation 
of the DO detector) are corrected according to the SME 
expectation of Eq. ([2|) . 

The SME contribution to the cross section has the 
general form /sME(i) = C^i,i?^(i)i?^(i)^"^ for the four 
model assumptions summarized in Table ID For each 
model, C^i, represents the constant coefficients we wish 
to determine, and A"^ refers to the appropriate linear 
combination of Ap and Ap . 



the different assumptions of Table U yields the following 
results: (i) Coefficients Ctt and Czz contribute only to 
the total cross section, and we do not attempt to extract 
these coefficients, (ii) Coefficients Ctx^, Cty, and Ctz 
combine with the small off-diagonal elements of matrices 
Ap and Ap , for which we expect poor sensitivity, (iii) 
Coefficients Cxz and Cyz couple to expressions that de- 
pend on sidereal time (differing by a phase of 7r/2). (iv) 
Coefficients Cxx and Cxy couple to time dependent ex- 
pressions with twice the sidereal frequency, and the two 
terms differ by a phase of 7r/4. 

Table [II] collects the resulting forms of the function 
/sme(0 foi" different assumptions. We refer to the "side- 
real phase" uJst as (p, where ojs is the inverse of the side- 
real day. The &-terms in these expressions depend on the 
colatitude of the detector, the orientation of the proton 
beam at the detector relative to geographic north, and 
the XX and ZZ elements of the combination of Ap and 
Ap that are appropriate to the particular assumption of 
the model. 



TABLE If: Forms for /sme((?!') used to extract SME 
coefficients. 



Condition 



/sme(0) 



Cxx 
Cxy ~ Cyx 
Cxz ~ Czx 
Cyz = CzY 



Cyy 2Cxx[ ^^^ cos 2<^ + 63 sin 24 



2CxY (^^^ sin 20 - 63 cos 2(? 
2Cxz (^4 cos 4> + bs sin 0) 
2Cyz (&4 sin 4> — b^ cos 0) 



Assuming that any LIV originates from just the top 
quark sector, we expect the background rate (principally 
VF-|-jets events) to be proportional only to the luminos- 
ity. To search for a signal varying with sidereal time, we 
sum the contributions to each of twelve Ni bins (corre- 
sponding to two sidereal hours each) for all data: 



TABLE L /sme(0 for different SME assumptions. 



Assumption 



/sME (t) 



(cu)^. = (cQ)^.i?^(f)i?^(i)(A^'° +^4^ 



(cq)^. =0 icu)^..Kit)mt){A 
d^,,R^^{t)Rf,{t)\Al' 



Q/Sn 



= 



c^^RmRmi^f + kAf) 



For each model, we estimate one possible component 
of Cpi/ at a time. We impose the requirements that each 
tensor C^i, is symmetric and traceless, choosing Cxx = 
—Cyy to satisfy the latter condition. We adopt the index 
ordering conventions n,v — {T, X, Y, Z} to refer to coor- 
dinates in the Sun-centered frame and a,(3 = {t^x^y,z} 
for coordinates in the DO frame. Evaluating Eq. ^ for 



N, « iVtot^ [1 + /s/sme(<^,:)] , (5) 

-'-int 

where iVtot is the total number of signal (ti) and back- 
ground (non-ii) events corresponding to the total inte- 
grated luminosity Ci^t, C.i is the integrated luminosity 
over the appropriate bin of sidereal phase (/)i, and ]§ is 
the average fraction of signal events in the data. 

We extract js from the data that was used previo usly 
to determine the tt cross section in ^-|-jets events [lO|. 
The ii cross section is measured in bins of jet multiplicity 
for the e-|-jets and /i-|-jets channels. The subset of events 
with at least four reconstructed jets that pass selection 
requirements contain a high fraction of ii events, provid- 
ing the best sensitivity to any time dependence in the 
tt event rate. We find /s(e-h>3-jets) = 0.78 ± 0.12 and 



/s(/x+>3-jets) = 0.76 ± 0.11. Because of this difference, 
we treat the electron and muon channels separately. 

To simplify fitting /sme('/>) to the data, we define a 
variable R for each bin: 



Ri — 



N,/Nu 

fs V ^^/^■ 



- 1 



int 



(6) 



Equation ^ is the luminosity-corrected sidereally- 
binned relative tt event rate, which can be compared 
directly to /sme(0)- In the absence of any significant 
sidereal time dependence, all the Ri values should be 
consistent with zero, while a sidereal time dependence 
would produce a sinusoidal variation in this rate. The 
amplitude for any sinusoidal dependence is given by the 
product of an SME coefficient and a mixture of contribu- 
tions from the rotation matrix and appropriate combina- 
tion of elements from Ap and Ap . This latter mixture 
also fixes the phase of the sinusoidal function in a fit to 
the data. 

The resulting distributions for _R as a function of side- 
real phase are shown in Fig.[l] separately for the electron 
and muon channels. The forms of /sme(0) are fitted to 
these two distributions to estimate the values of the SME 
coefficients for the assumptions summarized in Tables |T] 
andim We apply a smaU correction of 1.2%-4.7% to each 
extracted value to account for biases introduced by the 
finite bin size. 
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FIG. 1: The dependence of i?, as defined in Eq. ([6]), on 

sidereal phase for (a) e-|->3-jets tt candidates, and (b) 

^+>3-jets tt candidates. 

While the dominant contribution to the uncertainty 
on the SME coefficients results from the limited size of 
our tt data sample, the estimated fraction of ti events in 
the data contributes an additional uncertainty. We treat 
this as a systematic uncertainty in this study. The back- 
ground from single top quark events can, in principle, 
exhibit SME effects. However, their relative contribu- 
tion to the ti sample is negligible {^1%). The orienta- 
tion and location of the detector, as well as the origin 
chosen for the time of events, also carry negligible uncer- 
tainties. Finally, any uncertainties in the values of the 
elements of Ap and Ap can potentially contribute a 



systematic uncertainty in this analysis, in ways similar 
to those discussed in the analysis of the ti cross section, 
as summarized below. 

The leading sources of systematic uncertainty in the 
kinematics of ti events arise from: (i) the jet energy 
scale, (ii) jet energy resolution, and (iii) jet identifica- 
tion. These can affect the distributions of momenta re- 
constructed in the detector, but, as the elements of Ap 
and Ap reflect only average values of the components 
of the momenta over the detector acceptance, such aver- 
ages are not very sensitive to small changes in kinematic 
parameters. The relative uncertainty of the contributing 
elements is negligible compared to the statistical uncer- 
tainty of the data and the systematic uncertainties on 
signal fractions fs- 

A periodic time dependence could potentially be in- 
troduced to the event rate through changes in event se- 
lection efficiency. We check this possibility by exam- 
ining the luminosity-corrected sidereally-binned relative 
event rates {R distributions) for the lepton-|-n-jets chan- 
nels, where n = 2,3. These bins of jet multiplicity con- 
tain relatively small contributions from ti events, with 
/s(^+2-jets) « 12% and /s(^+3-jets) « 45%. We ex- 
tract the amplitudes for any time dependent oscillations, 
corresponding to the parameterizations used for the co- 
efficients in Tables IIIIfTVl in each of the four cross-check 
channels (^-|-n-jets, where ^ = e,/i and n — 2,3). For 
each assumption, the ensemble of fits is consistent with 
no time dependence at levels of probability in the range 
6%-38%. We therefore conclude that these cross checks 
give no indication of a sidereal time-dependent efficiency. 

Finally, it should be noted that any residual non- 
sidereal time dependence is suppressed greatly by folding 
the data into twelve bins of sidereal phase, as the mag- 
nitude of any residual contribution following this folding 
depends inversely on the difference in the period of the 
time-dependent efficiency and the sidereal period. Most 
problematic would be an unexpected time-dependent ef- 
ficiency with a period close to that of a sidereal day. 
The worst realistic case would be a contribution to de- 
tector efficiency that has a period of 24 solar hours. How- 
ever, because the data taking spans approximately seven 
years, any contributions from such an effect would be 
suppressed by about a factor of 10. To affect our conclu- 
sions, we would have had to experience a highly unlikely 
periodic dependence of the efficiency of approximately 
75% over 24 hours. No periodic effects of this magnitude 
have ever been observed in the detection efficiencies for 
objects considered in this analysis. 

Because the SME contribution to the matrix element is 
independent of lepton flavor, we perform a simultaneous 
fit to both the e-|->3-jets and /i-|->3-jets data to obtain 
the final results. The extracted SME coefficients are all 
consistent with no time dependence, and we therefore 
find no evidence for violation of Lorentz invariance in 
the ti system. 



We define the observed limits (95% C.L. intervals) for 
each SME coefficient as the extracted value ±2 stan- 
dard deviations. Because the magnitude of the 95% con- 
fidence bounds on elements of the linear combination 



^y.v 



(cq);ji/ + icu)fj,L' for the assumption of d 



liv 





are larger than 1, we cannot place meaningful limits on 
these combinations of SME coefficients in this analysis. 
The remaining limits are presented in Tables IIIII^Vl 



TABLE III: Limits on SME coefficients at the 95% 
C.L., assuming {cu)fj,i^ = 0. 

Coefficient Value ± Stat. ± Sys. 95% C.L. Interval 



{cq)xX33 


-0.12±0.11 


±0.02 


[-0.34, +0.11] 


iCQ)YY33 


0.12±0.11 


±0.02 


[-0.11, +0.34] 


{cq)xY33 


-0.04 ±0.11 


±0.01 


[-0.26, +0.18] 


{cq)xZ33 


0.15 ±0.08 


±0.02 


[-0.01, +0.31] 


(cq)yZ33 


-0.03 ±0.08 


±0.01 


[-0.19, +0.12] 



TABLE IV: Limits on SME coefficients at the 95% 
C.L., assuming {cq)^i, = 0. 



Coefficient 


Value ± Stat. ± Sys. 


95% C.L. Interval 


{cu)xX33 


0.10 ±0.09 ±0.02 


[-0.08, +0.27] 


{cu)yY33 


-0.10±0.09 ±0.02 


[-0.27, +0.08] 


(cu)xY33 


0.04 ±0.09 ±0.01 


[-0.14, +0.22] 


(cu)xZ33 


-0.14 ±0.07 ±0.02 


[-0.28, +0.01] 


{cu)yZ33 


0.01 ±0.07 ±< 0.01 


[-0.13, +0.14] 



TABLE V: Limits on SME coefficients at the 95% C.L. 



assuming c^i^ = 



Coefficient 


Value ± Stat. ± Sys. 


95% C.L. Interval 


dxx 


-0.11±0.10 ±0.02 


[-0.31, +0.09] 


dYY 


0.11 ±0.10 ±0.02 


[-0.09, +0.31] 


dxY 


-0.04 ±0.10 ±0.01 


[-0.24, +0.16] 


dxz 


0.14 ±0.07 ±0.02 


[-0.01, +0.29] 


dYZ 


-0.02 ±0.07 ±<0.01 


[-0.16, +0.13] 



In the SME, different particles can have distinct 
Lorentz-violating properties, so it is of interest to test 
all species. Most constraints on LIV are for particles of 
the first and second generations, with a few limits on 
SME coefficients for the third generation. The only sec- 
tor for which no constraints on Lorentz violation exist to 
date is the top quark [31 . The limits on the (cq)^i/33 and 
(c(7)^i/33 coefficients determined in this study represent 
the first constraints on LIV in the top quark sector, and 



the first such constraints on any free quark. 
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